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ABSTRACT 

We examine the physical properties of 173 cluster members in CL 1358+62 (z = 0.3283) from HST 
WFPC2 imaging taken in the F606W (^^rest-frame B) and F814W (~rest-frame V) filters over a 2.2 x 2.2 
Mpc^ field {Hq = 100 km s^^ Mpc~^, VIm = 0.3, f^A = 0.7). Structural parameters are measured for each 
galaxy by fitting a PSF-convolved, two component model to their 2D surface brightness distribution. We 
examine bulge+disk models using three different bulge profiles (de Vaucouleurs, Sersic, and exponential), 
and rigorously test the robustness of our results by analyzing several thousand artificial galaxies in the 
same manner as the cluster data. The measured physical properties from the best-fit profile of the cluster 
galaxies are combined with ground-based spectroscopy to test for correlations between morphological 
characteristics, current star formation, total galaxy colors, and cluster substructure. We find that: (1) 
Bulge-to-total ratio [{B /T)dev] and Hubble type (—5 < T < 8) are strongly correlated (99% confidence), 
but the scatter is large and early- type spirals are not reliably distinguished from ellipticals and SO's based 
on {B/T)deV- (2) From comparison of their physical properties, the low luminosity (—17.3 > MBz — 
51og/i > —19.3) ellipticals in our sample are likely to be face-on SO galaxies. (3) High galaxy asymmetry 
and strong [OIIJA3727 emission are strongly correlated for disk-dominated members [{B/T)deV < 0.4]. 
(4) There exists a small population (~ 5%) of bulge-dominated members whose significant [OIIJA3727 
emission (< — 5A) suggest they harbor active galactic nuclei. (5) At these redshifts, determining the 
correct Sersic index n can be highly unreliable. 



Subject headings: galaxies: clusters: individual: (CL 1358+62) 
galaxies: structure - galaxies: evolution 
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1. INTRODUCTION 

How galaxies evolve as a function of their environment 
remains a fundamental question in astronomy. In particu- 
lar, the environment of rich clusters provides a unique lab- 
oratory for studying how the physical properties of galaxies 
are related to varying local density, interactions with other 
galaxies, and exposure to a hot intracluster medium. One 
deeply ingrained and useful approach to studying galaxy 
evolution is characterizing morphological properties by vi- 
sually typing galaxies using the Hubble system (Hubble 
1926). As we move into the era of galaxy surveys with 
thousands, if not millions, of objects, however, there is a 
need for a quantitative, uniform, and reproducible method 
of cataloging these visual properties. 

A promising solution to this challenge is to measure 
physical properties, e.g. bulge and disk scale lengths, 
bulge-to-total ratios, disk inclination, bulge ellipticity, 
etc., directly from 2D surface brightness distributions. 



Here, parametric models are convolved with the point 
spread function and then compared to the galaxy's im- 
age to determine the best fit model (Schade et al. 1995; 
de Jong 1996; Marleau & Simard 1998; Simard et al. 1999; 
Tran et al. 2001; Trujillo et al. 2001; Simard et al. 2002; La 
Barbera et al. 2002; MacArthur et al. 2003). Not only are 
2D fits a natural way to weight the information contained 
in each pixel, they are crucial for detecting the presence of 
bars, HII regions, tidal tails, arms, shells, etc., i.e struc- 
tures not easily identified in traditional ID profiles. By 
keeping the 2D information, structural parameters for a 
large sample of galaxies can be derived and then com- 
pared to the galaxies' spectral properties and environmen- 
tal conditions. A lingering concern, however, is assess- 
ing how meaningful the structural parameters are. In the 
well- accepted framework of Hubble classifications, how can 
bulge-to-total ratios, scale lengths, galaxy asymmetry, etc. 
aid our understanding of galaxy evolution? 
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In Bubble's original system (Hubble 1926), one of the 
main qualitative classification criteria is bulge-to-disk ra- 
tio. As bulge+disk models measure a basic extension of 
this parameter (bulge-to-total fraction; B/T), it should 
be possible to correlate bulge fraction with Hubble type. 
The bulge fraction varies, however, depending on the type 
of bulge profile used, i.e. classical de Vaucouleurs r^/^ law 
(de Vaucouleurs 1948), the generalization by Sersic to r^/" 
(Sersic 1968), or pure exponential (de Jong 1996). A fun- 
damental question is which profile correlates best to Hub- 
ble classifications for galaxies spanning the range of mor- 
phological types? 

Once the appropriate profile is determined, the mea- 
sured structural parameters provide a wealth of informa- 
tion with which to study galaxies. In addition to confirm- 
ing the morphology-density relation in intermediate red- 
shift clusters (Dressier et al. 1997), recent studies find a 
population of late- type/disk-dominated members that lack 
strong emission (Poggianti ct al. 1999; Balogh et al. 2002), 
i.e. their star formation has been quenched. This result 
suggests that spectral and morphological evolution in the 
cluster environment may be decoupled and raises the ques- 
tion of how strongly ongoing star formation is correlated 
with physical characteristics. By comparing spectral fea- 
tures to quantitative measures of galaxy morphology from 
our extensive sample of CL 1358-1-62 {z = 0.3283; Fisher 
et al. 1998, hereafter F98), we can address this issue as 
well as test for correlations between morphology and clus- 
ter substructure {e.g. Hutchings et al. 2002). 

A direct test of the link between star formation and mor- 
phology is to compare [Oil] emission to residuals from the 
galaxy fits. Couch et al. (2001) find that spiral structure 
still can be observed in z ~ 0.3 cluster galaxies long after 
star formation has ceased, i.e. signatures of past events 
may be long-lived. If 1) morphological disturbances asso- 
ciated with star formation continue to exist long after the 
star forming episode ends, and 2) late-types are "stran- 
gled" when they enter the cluster environment (Balogh 
et al. 1998, R > 2h~^ Mpc), then we should observe a pop- 
ulation of predominantly late-type galaxies in CL 1358-1-62 
with high asymmetry but no detectable [OIIJA3727 emis- 
sion. By comparing [OIIJA3727 to bulge-to-total ratio and 
galaxy asymmetry, we can isolate these galaxies. 

Another interesting question that can be addressed us- 
ing 2D fits is whether low luminosity ellipticals and SO's 
are essentially the same type of galaxy. It has been argued 
that Mb > —22 ellipticals and SO's have different forma- 
tion/merger histories than the most luminous early- types. 
.Jorgcnscn & Franx (1994) suggest that these fainter el- 
lipticals and SO's actually form a common parent galaxy 
population, and that viewing angle plays a strong part in 
how these galaxies are typed (Rix & White 1990; Kor- 
mendy & Bender 1996, see also). Their similarity in col- 
ors support this argument (Sandage & Visvanathan 1978). 
As the isophotal shapes of early-type galaxies hold clues 
to their formation history (Nieto & Bender 1989; Naab 
et al. 1999), determining if bulge-dominated systems must 
be separated into two classes, i.e. "boxy" bright ellipti- 
cals and "disky" fainter ones (Kormendy & Bender 1996), 
is integral to understanding galaxy evolution across the 
range of luminosity (~mass). To address this question, 
we compare the bulge fractions, disk inclinations, bulge 



ellipticities, and velocity dispersions of visually classified 
ellipticals and SO's in CL 1358-^62 (Fabricant et al. 2000, 
hereafter FFvDOO). 

In this paper, we characterize the physical properties 
of 173 confirmed members in CL 1358+62 (F98) using 
GIM2D (Marleau & Simard 1998; Simard et al. 1999; Tran 
et al. 2001; Simard et al. 2002), an automated program 
that fits a PSF-convolved 2D surface brightness model to 
an image of each galaxy and searches space for the 
best fit. The images are from 8' x 8' HST WFPC2 mo- 
saics taken in the F606W and F814W filters (~ rest-frame 
B and V). The members span the range of Hubble type 
(—5 < T < 8) and have apparent F814W magnitudes be- 
tween 17.4 (BCG) and 23.2 (-20.8 > MB^ - 5 log /i > 
— 15.9). By fitting three different bidge+disk profiles to 
the sample, we determine the profile that is correlated best 
with the published visual classifications. We test for cor- 
relations between structural parameters, current star for- 
mation as defined by [OIIJA3727 emission, internal velocity 
dispersions, and local cluster structure. We also investi- 
gate the relationship between elliptical and SO galaxies, 
and examine the properties of E-|-A galaxies. 

To determine how robust our results are, we simulate 
thousands of artificial de Vaucouleurs bulge-|-exponential 
disk galaxies and fit surface brightness profiles to them 
in the same manner as with the cluster sample; note we 
assume galaxies can be well-described by a two compo- 
nent profile. The artificial galaxy catalog covers the same 
range in structural parameter space as the cluster mem- 
bers, and from these simulations we quantify systematic 
and random errors associated with bulge fraction, half- 
light radius, bulge & disk scale lengths, disk inclination, 
and bulge ellipticity. We also test our ability to measure 
Sersic and exponential bulge components by fitting artifi- 
cial galaxies with r^/" (0.2 < n < 6) bulge-l-exponential 
disk and double exponential profiles. 

The sections of this paper are organized as follows: The 
data and fitting technique are described in §2. In §3, we 
present the models and structural parameters determined 
by our surface brightness fitting program for the cluster 
sample. Here we employ "sanity checks" by comparing 
the measured bulge fraction to visual morphologies, lo- 
cal cluster density, and total galaxy colors. Correlations 
between bulge fraction, total residual, galaxy asymmetry, 
and current star formation also are examined in this sec- 
tion. We discuss the physical properties of elliptical, SO, 
and E-f-A galaxies in §4. In §5, we test the robustness of 
our measurements by utilizing results from thousands of 
artificial galaxies fits. Our conclusions are summarized in 
§6. Here we use Hq = lOO/i km s~^ Mpc~^, = 0.3, 
and = 0.7. 

2. DATA AND ANALYSIS 

2.1. Imaging 

We use HST WFPC2 imaging of CL 1358+62 that is 
part of an extensive project to study the evolution of in- 
termediate redshift cluster galaxies (van Dokkum et al. 
1998a,b; Tran ct al. 1999; van Dokkum et al. 1999, 2000; 
Kelson et al. 2000a,b,c, 2001; van Dokkum & Franx 2001). 
The 8' X 8' mosiac {R l.lh~^ Mpc) is centered on the 
brightest cluster galaxy (BCG) and is comprised of 12 sep- 
arate pointings in the F606W and F814W filters. Three 
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exposures, each 1200 seconds, were taken in each filter for 
eacli pointing. The two chosen filters, F606W and F814W, 
are similar to rest-frame B and V for CL 1358+62. The 
images were reduced at the Space Telescope Science Insti- 
tute (STScI) with the usual pipeline procedure. Details 
of the complex sky subtraction and cosmic ray removal 
are explained fully in van Dokkum et al. (1998b, here- 
after vD98). Following vD98, we transform from observed 
F606W and F814W magnitudes to rest-frame B and V 
using 

= F8UW + imi{F606W - F8UW) + 0.524 (1) 
= FSIAW + 0.204(F606T4^ - F8UW) + 0.652 (2) 

with a distance modulus of 40.03; we also correct for pas- 
sive evolution (van Dokkum et al. 1998a). 

2.2. Spectroscopy 

The large spectroscopic sample obtained for the 
CL 1358+62 field was taken at the Multiple Mirror Tele- 
scope and the William Herschel Telescope. A detailed de- 
scription of the reduction and analysis of the spectroscopic 
sample are described in Fabricant et al. (1991) and F98. 

Twenty slit-masks were designed to cover a 10' x 11' 
field; targets were chosen based on their R magnitude 
from imaging taken at the Whipple Observatory 1.2m 
telescope. The spectral resolution was ~ ISA at the 
WHT and - 20 A at the MMT. Of the 387 redshifts mea- 
sured, 232 are cluster members. Defining galaxies with 
0.31461 < z< 0.34201 as cluster members, CL 1358+62's 
mean rcdshift and velocity dispersion are 0.3283 ± 0.0003 
and 1027^11 km s^^ respectively (F98). 

From the rcdshift catalog, 192 cluster members fall on 
the 8' X 8' HST WFPC2 mosaic. The completeness within 
the HST mosaic is > 90% to i? = 21 (mg^ ~ 20.2), de- 
creasing to ^ 30% at i? = 22. The incompleteness is 
not due to an inability to measure redshifts at these mag- 
nitudes but by the limited number of galaxies observed. 
The redshift success rate is not strongly color dependent 
for i? < 23.5 (F98). 

Of the 192 galaxies that arc on the HST mosaic, we fit 
galaxy models only to 173 as 19 fall too close to a chip edge 
for proper analysis of the surface brightness distribution. 
The 173 galaxies range in absolute B^ magnitude from 
-20.8+51og/i mags (H375; BCG) to -15.9+51og/i mags 
(11829). We include the published [OII]A3727 A fine 
strength measurements (F98) and internal velocity disper- 
sions for a subset (Kelson et al. 2000b) in our analysis. 

2.3. Structural Measurements 

2.3.1. The Surface Brightness Models 

We use the GIM2D package (Marleau & Simard 1998; 
Simard et al. 1999; Tran et al. 2001; Simard et al. 2002) to 
find the best-fit PSF-convolved, 2D bulge+disk models to 
the surface brightness profiles of the cluster members. The 
program has a maximum of 12 fitting parameters: the flux 
{Ptotai) in the model integrated to r = oo; the bulge/total 
luminosity B /T = Ftuige/ Ftotai] the semi-major axis effec- 
tive radius of the bulge Tg; the bulge cUipticity e = 1 — b/a 
where a and b are the bulge semi-major and semi- minor 
axes respectively; the bulge position angle (j)b; the semi- 
major axis exponential disk scale length r^; the inclination 



of the disk i {i = for face-on); the disk position angle 
(j)d; the subpixel dx and dy offsets of the galaxy's cen- 
ter; the residual background level db; and the Sersic index 
n of the bulge. Both cf)), and (j)d are measured clockwise 
from the positive y-axis of the image. The best-fit param- 
eters and their confidence intervals are determined using 
the Metropolis algorithm (Metropolis et al. 1953; Saha & 
Williams 1994) which uses the xl test to determine the re- 
gion of maximum likelihood in the multi-parameter space. 
The bulge profile is defined as 

E(r) =Se ea;p{-fc[(r/re)i/"-l]} (3) 

where S(r) is the surface brightness at r along the semi- 
major axis, and Eg is the effective surface brightness. This 
bulge profile is also known as the Sersic profile (Sersic 
1968). The parameter k is equal to (1.9992n - 0.3271), 
a value that defines re to be the projected radius enclos- 
ing half of the light in the bulge component (Capaccioli 
1989). The classical de Vaucouleurs profile is a special 
case of Equation 3 with n ~ 4. 
The disk profile is defined as 

E(r) = So exp{-r/rd) (4) 

where Eq is the (face-on) central surface brightness. We 
note (as do Simard et al. 1999; Simard et al. 2002) that 
the bulge/disk nomenclature adopted here to describe our 
surface brightness models may not reflect the internal kine- 
matics of its components. A "bulge" may not be a cen- 
tralized, dynamically hot spheroid but could be a central 
star burst. Similarly, a "disk" may not necessarily be a 
cold, co-rotating population. For example, dynamically 
"hot" systems such as faint dwarf ellipticals are best flt by 
exponential disks (Binggeli & Cameron 1991; Ryden et al. 

1999) . 

Before a 2D model to the surface brightness of a given 
galaxy can be fitted, the galaxy's isophotal area and an 
appropriate point spread function (PSF) must be deter- 
mined. To define the isophotal area, we use the galaxy 
photometry package SExtractor V2.0 (Berlin & Arnouts 
1996) with a detection threshold of /i8i4 — 24.4 mags/D" 
(equivalent to ~ la of the sky noise) and a minimum de- 
blending contrast parameter of 0.01. As we fit galaxies 
in two filters (F606W and F814W), we use the isophotal 
area defined in the (redder) F814W image. To generate a 
PSF for each galaxy in each filter, we use Tiny Tim V4.4 
(Krist 1993). As the PSF changes across each WF chip, 
a PSF model is generated every 50 pixels and the nearest 
one to the galaxy is chosen for the GIM2D analysis. Both 
the PSF and galaxy models arc subsamplcd by a factor of 
five as WF image are undersampled and the pixelization 
can affect the shape of small galaxies such as those in our 
cluster sample. The centering of the galaxy is also im- 
proved by subsampling the data; this point can be critical 
for determination of the galaxy residuals (Conselice et al. 

2000) . 

By fitting models to the surface brightness distribution 

of these galaxies, we measure the structural properties n, 
B/T, Ve, Td, 4>bi (pd, i, and the half-light radius ri/2- The 
semi-major axis half-light radius is computed by integrat- 
ing the sum of Equations 3 and 4 to r = oo; note that as 
the models include bulge ellipticity and disk inclination, 
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circular symmetry is not assumed. The galaxy's (.x, y) cen- 
ter is also determined from the best fit model. 

The asymmetric image residual flux is quantified by the 
asymmetry index Ra (Schade et al. 1995), defined as 

Ra = {RA)raw — {RA)bkg 

= ii£ ^ (5) 

where Rij is the flux at (i, j) in the residual image, R}j'^ is 
the flux in the residual image rotated by 180°, and is the 
flux in the original image. Following Marleau & Simard 
(1998), Ra is measured within r = Ir^j^- The second term 
[RA)hkg in Equation 5 is a statistical correction for back- 
ground noise fluctuations. Since (RA)raw involves taking 
absolute values of pixel fluxes, it will yield a positive signal 
even in the sole presence of noise. The background cor- 
rection, (-Ryi)fcfeg, is computed over pixels flagged as back- 
ground pixels by SExtractor. The iJ^'s are background 
pixel values in the residual image, and the B\^^ 's are back- 
ground pixel values in the residual image rotated by 180°. 
The background correction is computed over a background 
pixel area equal to the pixel area over which (RA)raw is 
computed. Given the statistical nature of {RA)bkg, there 
will be cases when a galaxy is faint enough compared to 
its background noise that Ra raay take on small negative 
values in exactly the way as the difference of two values of 
{RA)bkg computed from different regions of the sky may 
be negative. 

In addition to Ra , wc measure the total residual fraction 
of light Rt by taking the pixels assigned to the galaxy by 
SExtractor and creating a mask that is applied to the orig- 
inal and model images (Tran et al. 2001). For a positive- 
definite residual fraction, the model is subtracted from the 
original, and the absolute value of the difference at each 
pixel over the isophotal area is summed. To account for 
the sky, the same number of sky pixels as galaxy pixels 
are summed in the same manner and subtracted from the 
total galaxy residual. We use the total residual fraction 
of light as a gauge of the model's goodness of fit but note 
that, for fainter galaxies, the error in Rt is dominated by 
the error in the sky flux. 

3. PHYSICAL PROPERTIES 

3.1. What Type of Profile? 

A fundamental question when using parametric surface 
brightness models to measure physical properties is what 
type of profile should be used. Depending on the pro- 
file, certain structural parameters {e.g. B/T, scale lengths, 
and bulge/disk colors) can change drastically in value (see 
Fig. 1). Recent work (de Jong 1996; Courteau et al. 1996; 
Andredakis 1998) combine two n = 1 exponential compo- 
nents to fit the surface brightness profile of late-type spi- 
rals while others (Caon et al. 1993; La Barbera et al. 2002) 
use a Scrsic profile to fit early- type galaxies. Another pop- 
ular profile is the classical de Vaucouleurs with exponential 
disk (Schade et al. 1995; Marleau & Simard 1998; Simard 
et al. 1999; Im et al. 2001; Tran et al. 2001). If conclusions 
are to be drawn from these quantitative studies of galaxy 



morphology, the most basic first step is to be consistent in 
the type of profile used. 

In this paper, we test three different profiles by fit- 
ting them to the cluster galaxies and comparing the mea- 
sured structural parameters to published Hubble types 
from FFvDOO; as we justify in §5, we use only the 155 
members in our sample brighter than msi4 = 21 [MB^ ~ 
—17.3 — 5-|-51og/i). The three profiles are: (1) de Vau- 
couleurs bulge with exponential disk; (2) Sersic bulge 
(0.2 < n < 6) with exponential disk; and (3) double expo- 
nential. Figure 1 shows B/T versus visual morphological 
type for the three profiles; here, visual types are repre- 
sented as E (-5), SO (-2), Sa (1), Sb (3), Sc (5), Sd (7), Sm 
(9), and Im (10); galaxies with high degree of asymmetry 
{Ra > 0.05, see §3.3.1) are noted with dots. 

Of the three profiles, both the de Vaucouleurs 
bulge-|-exponential disk and Sersic bulge-|-exponential disk 
correspond well to the visual types. The correlation be- 
tween measured B/T and visual type for both profiles 
is 99% (> 2a) significant with the Spearman rank test 
(Press et al. 1992). The trend between de Vaucouleurs 
bulge-l-exponential disk and type is not surprising because: 
1) bright ellipticals and the bulges of early- type spirals are 
well-fit by a de Vaucouleurs profile (Andredakis et al. 1995; 
Andredakis 1998); and 2) cluster galaxies tend to be bulge- 
dominated systems. We note, however, that there is large 
scatter in the correlation and that B/T cannot reliably 
distinguish early- type spirals from ellipticals and SO's. 

While the correlation between visual type and B/T for 
the Sersic bulge-|-exponential disk fits is strong (see Fig. 1), 
we do not recover the observed correlation between bulge 
power n and Hubble type found locally (Caon et al. 1993; 
Graham & Prieto 1999, see Fig. 2). This combined with 
results from the simulations discussed in §6 suggest that a 
Sersic bulge-|-exponential disk is not an appropriate model 
for this sample. Our sample does not contain enough spi- 
rals to test for the trend between n and early versus late- 
type spirals found by Graham (2001). 

There is no visible correlation between B/T and type 
for the double exponential profile. The mismatch between 
profile and galaxy type results in an alarming number of 
cluster galaxies (75%) being classified as disk-dominated 
sytems {B/T < 0.4). The lack of correlation is not sur- 
prising, however, given the inherent ambiguity between 
"bulge" and "disk" when describing the two components 
with the same functional form. Also, this profile is appro- 
priate mainly for late-type spirals (de Jong 1996; Courteau 
et al. 1996). These points convince us that the double ex- 
ponential profile is a poor choice for measuring structural 
parameters across the range of galaxy types. 

Based on these tests, we proceed to use the de Vau- 
couleurs bulge with exponential disk profile as our canon- 
ical fitting model for the sake of continuity across the full 
range of morphological types. 

3.2. 2D Surface Brightness Models: Sanity Checks 

Before drawing any conclusions from the surface bright- 
ness models, we employ "sanity checks" to ensure our 
method corresponds to observed, well-established cluster 
galaxy properties. These tests are: 1) check for consis- 
tency of B/T measurement between filters; 2) establish a 
correlation between measured B/T and visual morpholog- 
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ical type; 3) determine if a morphology-radius (~ density) 
relation using B/T exists; and 4) check that typing cluster 
members hy B/T results in separation of galaxy types in 
the color-magnitude diagram. 

3.2.1. Galaxy Models: de Vaucouleurs Bulge with 
Exponential Disk 

We fit a do Vaucouleurs bulge with exponential disk 
model to the 173 cluster members in our sample. The 
galaxies are fit in each filter separately to test the robust- 
ness of the structural parameter measurements. Figure 3 
shows 168 of the 173 galaxies in our cluster sample. Each 
set of three thumbnails shows the galaxy, its best-fit de 
Vaucouleurs with exponential disk model, and its residual 
image created by subtracting the model from the original. 
Included in the images arc reference numbers (upper left) , 
apparent F814W magnitudes (bottom left), bulge/total lu- 
minosity {B/T; bottom middle), and the asymmetry pa- 
rameter and total fraction of residual light [RatRt] bot- 
tom right). The pixel area in each thumbnail is 15 times 
the galaxy's isophotal area as defined by SExtractor; this 
is usually 10 — 15" on a side. The cluster members range in 
absolute B^ magnitude from -20.8-F51og/i (BCG; H375) 
to — 15.9-|-51og/i (H1829). Table 1 lists these parameters as 
well as Bz, {B — V)z, Hubble type, half-light radius, and 
bulge/disk scale lengths. In the following discussion, we 
only include the 155 members with msu < 21 as, based 
on results from extensive simulations in §5, our measured 
structural parameters arc not significantly biased for these. 

For simplicity, we focus on structural parameters mea- 
sured from the individual F814W (~rest-frame V) images 
as opposed to including values from the F606W (^ rest- 
frame B) images. Comparing structural parameters mea- 
sured in the two filters finds little difference in their general 
distributions. We emphasize that while structural parame- 
ter values for individual galaxies may vary between the two 
filters, these differences are 1) expected and understand- 
able and 2) do not impact the global conclusions drawn 
from our analysis. 

3.2.2. B/T vs. Visual Morphology 

In Fig. 4, we show the bulge distributions measured in 
F814W. Not surprisingly, the majority of cluster mem- 
bers (~ 70%) are bulge-dominated {B/T > 0.4) systems. 
This result does not depend strongly on filter as, using the 
Kolmolgorov-Smirnov test (Press et al. 1992), we find the 
F814W and F606W distributions to be indistinguishable. 

With the de Vaucouleurs bulge+exponential disk pro- 
file, we find a strong trend (99% confidence) between B/T 
and visual morphological type (Fig. 1, top panel) with the 
Spearman rank test (Press et al. 1992), albeit with large 
scatter. We also see in Fig. 1 that of the high asymmetry, 
visually typed galaxic;s, all have B/T < 0.7; this is ex- 
pected because disks are more likely to harbor asymmet- 
ric features such as HII regions. The correlations between 
B/T, asymmetry, and visual type are highly encouraging 
but we note that, in this sample, it is difficult to sepa- 
rate early-type spirals from ellipticals and SO's using B/T 
alone. 

3.2.3. B/T (Morphology) vs. Local Density 



F98 and vD98 have shown that galaxy spectral type 
and visual morphological type vary as a function of local 
density in CL 1358-1-62. In Fig. 5, we show that fitting a 
de Vaucouleurs bulge with exponential disk profile to the 
same sample of galaxies also recovers the same results, as 
expected since we have demonstrated the strong correla- 
tion between B/T and visual type. Splitting the sample 
into bulge and disk-dominated systems, Fig. 5 shows the 
increase in the fraction of disk-dominated systems with 
decreasing local galaxy density. Here, the local galaxy 
density is defined as S = ll/Trrfg, and rio is the distance 
to the farthest of the ten nearest confirmed members to 
the object (Dressier 1980). 

An interesting question is whether the distribution of 
bulge-dominated galaxies is related to cluster substruc- 
ture. Using the Dressler-Shectman test (Dressier & Sheet- 
man 1988), F98 found the degree of substructure to be sig- 
nificant with 96% confidence (> 2a) where the significance 
was calibrated with 10^ Monte Carlo realizations. Fig. 6 
shows the spatial distribution of m8i4 < 21 mags members 
as a function of B/T and 6, the substructure statistic. As 
5 quantifies how much the local mean redshift and velocity 
dispersion deviate from the cluster's global values, groups 
of large circles indicate a significant degree of substructure. 
Moving from bulge to disk-dominated systems, we find 
that bulge- dominated systems are clustered more strongly 
than disk-dominated ones. The spatial distribution of the 
most bulge-dominated galaxies {B/T > 0.75) is different 
from that of the most disk-dominated {B/T < 0.25) ones 
with 97% confidence (> 2a) using the two dimensional K-S 
test. The most significant subclump is a foreground group 
composed mainly of early- types {B/T > 0.4), including 
the BCG. 

As a further refinement, we separate the cluster sam- 
ple into B/T > 0.5 and B/T < 0.5 members and apply 

the Dresslcr-Schectman test; again, the substructure sign- 
ficance is calibrated with 10"^ Monte Carlo realizations. 
The bulge-dominated galaxies (98) show substructure with 
98% confidence (> 2a) while the disk-dominated galaxies 
(75) lack significant substructure. In addition, the cluster 
velocity dispersion determined from bulge-dominated sys- 
tems is lower than that from disk-dominated ones: 950±70 
km s~^ compared to 1220 ± 100 km s~^ . 

3.2.4. Color- Magnitude Relation 

Following vD98, we examine how galaxy colors change 
as a function of magnitude and galaxy type (as defined 
by bulge fraction B/T). In Fig. 7, we show the color 
magnitude diagram in four B/T bins; we use B^ magni- 
tudes, {B — V) colors, CM relation (determined from E-SO 
members) from vD98. As B/T decreases, the scatter and 
average deviation from the CM relation increases. Also, 
the most disk-dominated systems {B/T < 0.25) are signif- 
icantly bluer and have the largest scatter about the CM 
relation. These results are consistent with the trends vD98 
found using Hubble classifications for the same galaxies. 

In Fig. 7, we see there exist cluster members with signif- 
icant disk components {B/T < 0.5) that are also red. It is 
possible these disk systems have had their star formation 
"strangled" Balogh et al. (1998) and/or they have man- 
aged to avoid major galaxy interactions that would have 
disrupted their disks, e.g. mergers. An extensive analysis 
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of the bulge/disk colors of these galaxies will be discussed 
in a future paper. 

3.3. Image Residuals of Cluster Members 

3.3.1. High Asymmetry Galaxies 

Following Tran et al. (2001) and Schade et al. (1995), 
we define galaxies with Ra > 0.05 as having a high de- 
gree of asymmetry. In the cluster sample of galaxies with 

msu < 21, ^ 10% have Ra > 0.05, and all of these high 
asymmetry galaxies have B/T < 0.7. Most asymmetries 
in this population are due to star-forming regions in disk- 
dominated systems {e.g. H178, H200, H234, and H396; see 
Fig. 8). 

CL 1358+62's fraction of high asymmetry galaxies is 
more similar to that observed in groups than in the 
field. Tran et al. (2001) found that about 11% of galax- 
ies nearby, X-ray luminous groups have high asymmetry^, 
with the fraction of bulge-dominated galaxies in groups at 
^ 50 — 60%. In comparison, Schade et al. (1995) found 
that 10/32 (30%) HST observed field galaxies within the 
redshift range 0.5 < z < 1.2 have high galaxy asymme- 
try. This is as expected, however, as the majority of field 
galaxies are disk-dominated, and disks are more sensitive 
to morphological disturbances. 

3.3.2. High Total Residual Galaxies 

After examining the residual images (Fig. 3), we define 
galaxies with Rt > 0.1 as having high total residual. Ap- 
proximately 20% of the members have high total residu- 
als with features that are over or under-subtracted by the 
model, e.g. bars (H525, see Fig. 8) or star-forming regions 
(H178, see Fig. 8). As there are 15 galaxies common to 
both populations, however, this is not surprising. 

We find that disk-dominated galaxies tend to be more 
disturbed than bulge-dominated ones: a third of these 
systems show high asymmetry and half have a high to- 
tal residual. Comparatively, the early-type galaxies have 
smoother galaxy profiles: only ^ 5% show high galaxy 
asymmetry and ~ 10% have a high total residual. Nor- 
malizing the residuals by the average surface brightness 
instead of the total flux finds the same result. The differ- 
ence in the high residual fraction of disk- dominated mem- 
bers compared to bulge-dominated ones supports findings 
that the mechanism(s) responsible for morphologically dis- 
rupting cluster members are: (1) more effective at disturb- 
ing disk-dominated galaxies than bulge-dominated ones 
and/or (2) the early- type galaxies are older and have been 
in the cluster environment longer, so any disruptions in 
these galaxies occurred long ago and have damped out via 
phase mixing. 

3.4. Correlations between Structural Parameters and 
[On]X3727 Emission 

To determine how morphological characteristics arc cor- 
related (or not) with ongoing star formation, we compare 
B/T, Ra, and Rt to [OII]A3727 measurements from F98. 
We assume [Oil] emission is correlated with star forma- 
tion and not active galactic nuclei (AGN) as we do not 
have the necessary data (e.g. Ha or X-ray) to distinguish 
between the two. Again, we discuss only the 155 galaxies 



in the cluster sample with msu < 21. In the following 
discussion, note that neither B/T nor Ra are correlated 
with luminosity, as determined with the Spearman rank 
test. 

3.4.1. Current Star Formation vs. B/T 

Figure 9 shows the distribution of the galaxies' bulge 
fractions versus their [Oil] equivalent widths; negative 
[Oil] EW values correspond to emission. Of the 26/155 
galaxies with significant [Oil] emission ([Oil] EW< -5A), 
~ 60% arc disk-dominated galaxies {B/T < 0.4). Even 
more compelling, 4 of the 7 galaxies with the highest star 
formation rates ([OII]< — 20A) have B/T = 0, i.e., com- 
pletely disk- dominated systems. 

Current star formation in cluster galaxies, as traced 
by [OII]A3727 emission, is not confined to only disk- 
dominated systems (see bottom panel. Fig. 9). Inter- 
estingly, ^ 40% of galaxies with significant [Oil] emis- 
sion are bulge-dominated {B/T > 0.4). These galaxies 
span the range of luminosity (~ mass), having —17.3 > 
MB^ - 5 log/i > -20.8, and include the BCG. We discuss 
this population in more detail in §4.2. 

Comparison of the [Oil] emission line strength to bulge 
fraction finds it to be correlated strongly with B/T. Ta- 
ble 2 shows the average [Oil] EW in four B/T bins. The 
average [Oil] emission is dramatically higher for galax- 
ies with B/T < 0.25, dropping to zero for B/T > 0.85 
galaxies. Like Balogh et al. (1998), we find that the most 
bulge-dominated members {B/T > 0.8) have [Oil] EQW 
consistent with zero. 

3.4.2. Current Star Formation vs. Galaxy Residuals 

The correlation between high galaxy asymmetry and 
strong [OII]A3727 emission for bulge and disk-dominated 
galaxies is shown in Fig. 10. Approximately 70% of the 
high asymmetry galaxies have significant [Oil] emission. 
For the disk-dominated members, there is a strong trend 
(99% confidence) of increasing [Oil] emission with in- 
creasing Ra] most of these galaxies have high asymmetry 
and/or high total residuals and are blue. We find that 
current star formation in the disk- dominated population 
is highly correlated with measureable morphological dis- 
turbances. 

Although 11 of the bulge-dominated galaxies do have 
strong [Oil] emission, only three of these galaxies have 
Ra > 0.05 and/or Rt > 0.1. Assuming the [Oil] emission 
is due to star formation, active star formation tends to be 
distributed in a smooth and uniform manner when it oc- 
curs in B/T > 0.4 members. Alternatively, these members 
may harbor active galactic nuclei (Martini et al. 2002) ; we 
discuss this possibility in greater detail in §4.2. Using the 
high total residual population instead of the high asym- 
metry one, we find essentially the same results. 

4. TOWARDS UNDERSTANDING DIFFERENT GALAXY 
POPULATIONS 

4.1. E's vs. SO's: A Matter of Viewing Angle? 

Jorgensen & Franx (1994) suggest that Mb > —22 ellip- 
ticals and SO's form the same galaxy class, such that ellip- 
ticals tend to be face-on members of this class while SO's 



Note, however, that our HST data are ~ 3x higher in physical resolution than the ground-based sample (0.44?i ^ kpc versus 1.31/1-1 kpc). 
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are more edge-on members. Our measurements of bulge el- 
lipticity and disk inclination enable us to test this hypoth- 
esis. If the ellipticity and inclination distributions of these 
two populations are completely different but they span 
the same range in other physical properties, this would 
be strong evidence for these ellipticals being face-on SO's. 
However, as the most luminous E/SO galaxies {e.g. BCGs) 
are hypothesized to have different formation histories (Ko- 
rmendy & Bender 1996; Naab et al. 1999), we restrict our 
analysis to 19 < < 21. 

Classically, ellipticals can be circular (EO; axis ratio is 
1) to cigar-shaped (E6; axis ratio is 0.4). Thus, visually 
typed elhpticals in CL 1358+62 (FFvDOO) should span 
the range in bulge ellipticity and, for systems with disks, 
inclination; the same is true for SO and even Sa galax- 
ies. To test this, we plot in Fig. 11 the bulge ellipticity 
and disk inclination versus measured B/T of E's and SO's. 
We include all visually typed E's and SO's in Fig. 11 for 
completeness but note that ellipticity and inclination mea- 
surements are robust only for B/T > 0.4 and B/T < 0.6 
systems respectively. 

Applying a K-S test to both the bulge ellipticity and 
disk inclination distributions of E's and SO's confirms they 
are different, yet their B/T distributions (for B/T > 0.4) 
are indistinguishable. Note also the conspicuous lack of 
round {Ell < 0.2) SO's, a population that should exist due 
to projection effects alone. To strengthen our argument, 
wc! compare the internal velocity dispersion'', half-light ra- 
dius, projected cluster radius, [OIIJA3727, magnitude, and 
{B - V)^ distributions of the E's and SO's. Using the K-S 
test, we find the distributions to be indistinguishable for 
all these parameters. 

The fact that visually typed E's and SO's share simi- 
lar properties with the exception of bulge ellipticity and 
disk inclination strongly suggests that these ellipticals are 
likely to be face-on SO's. In addition, the range in B/T 
for members classified as ellipticals (Fig. 11, top panel) 
shows that these are not pure r^/^ galaxies. This supports 
results from recent observations and modeling (Rix ct al. 
1999; Burkert & Naab 2001) that a disk component exists 
in most ellipticals. 

4.2. AGN's in Cluster Ellipticals? 

We find ^ 5% of the total cluster population are 
bulge-dominated galaxies {B/T > 0.4) with significant 
[OII]A3727 emission (< -5A) that, unlike their disk- 
dominated counterparts, display little morphological dis- 
ruption {Ra < 0.05). A possible explanation for these 
unusual systems is that they are similar to the blue field 
ellipticals studied by Menanteau et al. (1999, 2001), i.e 
bulge- dominated systems still forming stars in their cores. 
In this cluster, however, the majority of these systems are 
as red as the absorption population (see Fig. 10) even 
though their average [OII]A3727 emission is significantly 
higher (-8.5 ± 3.9 A vs. 0.70 ± 1.6A). Their bulge compo- 
nents are equally as red, i.e. these objects do not have blue 
cores and it is highly unlikely that their [Oil] emission is 
due to a central star burst. 

Our sample finds that the undisturbed, high [Oil] emis- 
sion, bulge-dominated members arc not low mass objects. 
Comparison of internal velocity dispersions measured for 

^ Prom Kelson et al. (2000b) we have a^igp for four ellipticals and 11 



half of this population (6/11) from Kelson et al. (2000b) 
shows (Tdisp ranges from 128 ± 5 km/s to 230 ± 5 with the 
BCG at 307 ± 7 km/s. While it could be argued that 
velocity dispersions are measured most easily for mas- 
sive galaxies, thus introducing a bias exists against low 
mass (~ (Jdisp) objects, the sample published by Kel- 
son et al. (2000b) includes robust velocity dispersions to 
(T 60km/s. 

These points suggest that the [Oil] emission in these 
red, bulge-dominated galaxies originates from AGN ac- 
tivity rather than star formation. In Abell 2104, Mar- 
tini et al. (2002) find a surprisingly high number of X-ray 
sources associated with red cluster galaxies; their lower 
limit on the AGN fraction of 5% is remarkably similar 
to our cluster fraction of low asymmetry, [Oil] emitting, 
bulge-dominated members. Unlike the AGN candidates 
in this cluster, however, the majority of Martini et al. 
(2002) 's AGN candidates lack significant optical emission. 
It may be that the number of AGN candidates identified 
via [Oil] emission only places a lower limit on the cluster 
AGN fraction. 

4.3. E+A Galaxies 

Following F98, we select post-starburst ("E+A") galax- 
ies as having no [OII]A3727 emission (> — sA) and strong 
Balmer absorption [(H/?+H(5+H7)/3 > 4 A). The nine 
spectroscopically confirmed E+A galaxies on the HST mo- 
saics span the range in bulge fraction (0 < {B/T) < 0.9) 
but these are not bright galaxies (—19.6 < MB^ — 51og/i < 
— 17.3). Virtually all of them are blue (see Fig. 7, also 
Bartholomew et al. 2001), and a third of them show sig- 
nificant asymmetry or total residual (see Fig. 10). A de- 
tailed analysis of the cluster E+A population is presented 
in Tran et al. (2003). 

5. SIMULATIONS: ROBUSTNESS OF SURFACE BRIGHTNESS 

MODELS 

To test the robustness of structural parameters deter- 
mined by fitting de Vaucouleurs bulge+exponential disk 
surface brightness models, we generate several thousand 
artificial galaxies and apply the same analysis. Note that 
we assume throughout this section that galaxies are de- 
scribed exactly by an r^/" bulge+exponential disk profile. 
In the following, we focus on errors associated with in- 
correct measurements of flux, B/T. bulge ellipticity, disk 
inclination, half-light radius, and the scale lengths of the 
two profile components. We also include short discussions 
on the reliability of fitting Sersic and exponential bulges. 

5.1. Determining de Vaucouleurs Bulges 

Artificial de Vaucouleurs bulge+exponential disk galax- 
ies are generated uniformly to cover the same range in 

flux, bulge/total luminosity, scale lengths (r^, r,.), bulge 
eccentricity, and disk inclination as the cluster sample in 
the F814W filter (3600s integration). For simplicity, the 
bulge and disk position angles are fixed at 45°. Once the 
smooth galaxy image is created, it is convolved with the 
appropriate TinyTim PSF, Poisson noise is added, and it 
is imbedded in a 30" x 30" WF2 image section that is 

I's; for the rest, we use derived velocity dispersions (Tran et al. 2003). 
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free of any detectable objects. At this point, the artificial 
galaxy catalog is analyzed in exactly the same manner as 
the observed sample. 

By comparing measured to input values of the structural 
parameters for the artificial galaxies, wc can map system- 
atic and random errors due to the galaxy's position in 
the multi-dimensional structural space. Systematic error 
is the median difference between the input and measured 
value for a parameter, and random error is the la width 
of the distribution associated with the median difference. 
Statistics drawn from the simulations are robust across 
the range of msu and {B/T)su as the artificial catalog is 
well-sampled in each of these bins. 

5.1.1. Total Magnitude 

The top panel of Fig. 12 shows the median fractional 
difference between input and measured magnitude (F814) 
over the range of B/T for the artificial catalog. For the 
brightest galaxies, the measured flux does not deviate sig- 
nificantly from the input value. At msu > 22 mags, 
however, we systematically underestimate the input fiux 
by ~ 0.3 mags and maximum random errors increase to 
Am8i4 > 0.5 across the range of B/T. With the exception 
of three galaxies (of 173), all of the CL 1358-1-62 sample 
are brighter than msu ~ 22 mags. 

The bottom panel of Fig. 12 emphasizes how faint galax- 
ies, and in particular low surface brightness ones, are ad- 
versely affected by decreasing signal-to- noise ratios. In 
each subpanel, wc underestimate the total flux as mea- 
sured ri/2 increases, i.e. as the average surface bright- 
ness decreases. Even at fairly bright magnitudes {msi4 < 
21 mags), flux loss can be significant for galaxies with 
ri/2> 5h~^ kpc. Fortunately, the median half-light ra- 
dius of the CL 1358-1-62 cluster sample is 2.4± 1.2/i~^ kpc 
so this should not affect our conclusions. 

5.1.2. Bulge/Total Luminosity 

Figure 13 illustrates how the measured bulge fraction 
[[B /T)out] differs from its input value as a function of 
measured rusu, B/T, and ri/2. Above our magnitude cut 
of msi4 = 21, the systematic differences between input 
and measured B/T are negligible and the random errors 
are < 0.2 (sec Fig. 14). These results lend confidence to 
measurement of the bulge/total luminosity for galaxies in 
our cluster sample. 

5.1.3. Half-Light Radius and Bulge/Disk Scale Lengths 

Figure 15 (top panel) shows the fractional difference be- 
tween input and measured values of ri/2 as a function of 
bulge fraction. Recovery of the half-light radius is robust 
for 771814 < 21 mags galaxies with systematic and random 
errors of less than 10% and 20% respectively. This is re- 
iterated in the top panel of Fig. 16. Here we show the 
average input and measured values as a function of 
input bulge fraction. Only for the most bulge-dominated 
systems [B/T > 0.7) are the differences systematically 
larger than la as due to the long wings of the de Vau- 
couleurs profile, flux at large radii can be lost in the sky 
background. 

Recovery of disk exponential scale length for B/T < 0.4 
galaxies is excellent. At m8i4 < 21 mags, systematic and 



random errors for disk-dominated objects {B/T < 0.4) 
are zero and < 15% respectively (middle panel. Fig. 15). 
In fact, the simulations show that measurements are 
robust up to B/T ~ 0.8 (middle panel. Fig. 16). 

In comparison, measurements of the bulge scale length 
are not as robust: re tends to be underestimated (~ 
10 — 20%) at all magnitudes (bottom panel. Figs. 15 & 
16). Random errors, even for the most bulge-dominated 
systems, also are larger than those associated with the 
disk scale length. While this result is not surprising when 
considering the shape of the de Vaucouleurs profile, it can 
introduce a bias that would be particularly detrimental to 
large, bulge-dominated, low luminosity systems. 

From these simulations, we determine that for the 
CL 1358-1-62 cluster sample: 1) Measurements of half- 
light radius across the range ol B/T are robust; 2) Disk 
scale lengths are reliable for galaxies with substantial disks 
{B/T < 0.6); and 3) The systematic uncertainties associ- 
ated with measuring bulge scale lengths requires that care 
be applied in their interpretation, especially for galaxies 
with a small bulge component {B/T < 0.4). 

5.1.4. Bulge Ellipticity & Disk Inclination 

To test the robustness of our results concerning ellipti- 
cals and SO's, we compare input and measured values of 
bulge ellipticity and disk inclination. Figure 17 shows the 
average differences in bulge ellipticity and disk inclinations 
for artificial galaxies as a function of measured B/T and 
apparent magnitude. At m8i4 < 21, both parameters are 
recovered well with median differences between input and 
measured values of approximately zero. The associated 
random errors in both ellipticity and cos i are < 0.1. As 
expected, random errors in the bulge ellipticity are largest 
for disk-dominated systems, and errors in cos i largest in 
bulge-dominated systems. 

5.2. Determining Sersic Bulges 

We find that fitting Sersic bulges with exponential disks 
to the cluster members results in a significant correlation 
between visual type and bulge fraction (see Fig. 1). By 
fitting Sersic profiles, however, we must test our ability to 
recover the true bulge power n. To do so. we create an ar- 
tificial galaxy catalog of 1500 galaxies in the F814W filter 
with bulge fraction > 0.5 and bulge power n between 0.2 
and 6; all other structural parameters are the same as for 
the de Vaucouleurs bulge-|-exponential disk artificial cata- 
log. As before, we discuss results only for galaxies with 
TO814 ^ 21 mags. 

5.2.1. Recovery of the Sersic Index n 

For most galaxies (70%), we are within 0.5 of the true 
Sersic bulge power n (Fig. 18; top panel). Approximately 
20% of the catalog, however, has An > 1 and, in some 
cases. An can be as much as ~ 4. Even more discon- 
certing is the comparison of input versus measured bulge 
fraction (Fig. 18; lower panel). The distribution is skewed 
heavily towards underestimating the true bulge fraction 
(median A{B/T) ~ 14%), and 23% of the galaxies are 
now considered disk-dominated systems even though none 
have input B/T < 0.5. 

These simulations emphasize that measuring reliable 
Sersic n and the corresponding bulge fraction requires 1) 
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a combination of excellent resolution and high signal-to- 
noise; 2) understanding how noise affects the measure- 
ments; and 3) rigorous testing of structural parameters 
derived for each galaxy. 

5.2.2. Fitting Sersic Bulges with r^/^ Profiles 

With the artificial catalog of Sersic bulge+exponential 

disk galaxies, we address the question of how serious the er- 
rors would be if we modeled them as de Vaucouleurs (r^/'*) 
bulges. Figure 19 compares the input Sersic B/T com- 
pared to that recovered using the de Vaucouleurs profile. 
The distribution is as skewed towards underestimating the 
true bulge fraction as in Fig. 18; this is due primarily to 
high B/T systems with bulge n ~ 1 being confused as 
disk-dominated systems (see §5.3). Even for galaxies with 
higher Sersic indices (n > 3), the average A{B/T) is still 
larger than if the true bulge profile is r^/*. 

However, we find that the half-light radii (^sizcs) are 
robust. As for bulge ellipticities and bulge scale lengths 
of bulge- dominated [{B/T)out > 0.4] galaxies, the sys- 
tematic errors are still negligible but the random errors 
are larger by ~ 2 and ~ 3 respectively. Despite the un- 
certainty in (B/T), certain structural parameters such as 
half-light radius as well as bulge ellipticity and scale length 
[for {B/T)ouT > 0.4] are reliable. 

5.3. Determining Exponential Bulges 

Although the correlation between visual type and bulge 
fraction as determined using a double exponential profile is 
poor (see Fig.l), we test our ability to recover this profile 
for completeness. Wc create an artificial galaxy catalog 
of 1500 galaxies with exponential bulges and disks in the 
F814W filter that spans the range of B/T. Again, we limit 
our discussion to galaxies with nisu < 21 mags. 

Figure 20 illustrates how the bulge/disk components 
often are reversed suc;h that bulge-dominated systems be- 
come disk-dominated ones and vice versa. This ambiguity 
manifests as the long tails in the two panels of Fig. 20. 
While our inability to recover the true double exponen- 
tial profile is not a compelling argument for not using it, 
these simulations demonstrate that a more sophisticated 
version of this profile is needed. For example, incorporat- 
ing a truncation radius (van der Kruit & Searle 1982; van 
dcr Knut 1987; de Grijs et al. 2001; Krcgel ct al. 2002) for 
the bulge component may alleviate the ambiguity. 

6. CONCLUSIONS 

Combining deep wide field HST imaging (8' x 8' field; 
2.2 X 2.2 Mpc^) with ground-based spectroscopy, we char- 
acterize the physical properties of 173 cluster members in 
CL 1358-h62 {z = 0.3283). Our study is unique due to 
the large number of confirmed cluster members, the dis- 
tance to the cluster, and the high physical resolution of the 
galaxies (0.27/i~^ kpc). By fitting two component surface 
profiles directly to each galaxy, we measure their structural 
parameters and compare these results to current star for- 
mation, total galaxy colors, and cluster substructure. In 
this study, we examine the viability of three different bulge 
profiles (de Vaucouleurs, Sersic, and exponential) paired 
with an exponential disk component. 

We stringently test the robustness of our results by fit- 
ting thousands of artificial galaxies in the same manner 



as the sample; note that we assume galaxies are well- 
described by an r^/" bulge-|-exponential disk profile. As 
the artificial galaxy catalog spans the same range in lumi- 
nosity, bulge fraction, half-light radius, and scale lengths 
as the cluster sample, we can quantify the systematic and 
random errors associated with the structural parameter 
measurements. From these simulations, we find that we 
tend to underestimate galaxy bulge fraction, especially for 
bulge-dominated objects. However, over the selected mag- 
nitude range (r«,8i4 < 21 mags; MB^ < — 17.3+51og/i) for 
the de Vaucouleurs bulge-|-exponential disk profile, mea- 
surements of B/T have systematic errors of < 15%. 
Our conclusions are as follows: 

1. In this sample, the correlation between {B/T)dev 
and Hubble type (-5 < T < 10) is strong (99% 
confidence), but the scatter is large. In particular, 
early-type spirals cannot be distinguished from el- 
lipticals and SO's based on their bulge fraction. 

2. Using quantitative morphological parameters, we 
find a predominance of bulge-dominated (early- 
type) galaxies and an established morphology- 
density relation, both characteristics of a well- 
developped cluster. 

3. The physical properties, i.e. bulge ellipticities, disk 
inclinations, B/T, half-light radii, and internal ve- 
locity dispersions, of the low luminosity (—17.3 > 
MB^ - 51og/i > -19.3) ellipticals and SO's in our 
sample suggest that these ellipticals are likely to be 
face-on SO galaxies. 

4. We find a population (~ 5%) of bulge-dominated, 
massive [a > 130 km ) members with signifi- 
cant [Oil] A3727 emission (< — 5A), morphologically 
smooth profiles, and colors as red as the absorp- 
tion galax;ies; the latter is true for both their total 
and bulge component. These physical characteris- 
tics suggest they harbor AGN. 

5. The most disk-dominated members {B/T < 0.25) 
have significantly higher average [OII]A3727 equiva- 
lent widths than members with B/T > 0.25. Disk- 
dominated galaxies also tend to have high degrees 
of structural asymmetry. 

6. Mimicking the sample with thousands of artificial 
galaxies is integral in understanding how robust cor- 
relations between physical and spectral properties 
are. In particular, measuring accurate Sersic in- 
dices at these redshifts requires significantly higher 
signal-to- noise and/or resolution than that of our 
WFPC2 imaging. 
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Fig. 1. — Comparison of bulgc-to-total fractions to Hubble types from FFvDOO for 155 typed cluster members {msi4 < 21, MBz ~ 
— 17.3+51og/i). The top panel shows B/T as defined using a de Vaucouleurs bulge+exponential disk, the middle a Sersic bulge (0.2 < n < 
6)+exponential disk, and the bottom a double exponential profile; all profiles are measured in the F814W filter (~rest-frame V). The small 
open circles represent E-SO's (—5 < T < —1), open triangles SO-a (0 < T < 1), open squares spirals and irregulars (2 < T < 15), and large 
open circles post-starburst galaxies; members with a high degree of galaxy asymmetry {Ra > 0.05) have a solid dot. The vertical line denotes 
the adopted break between bulge (B/T > 0.4) and disk (B/T < 0.4) dominated galaxies (Tran et al. 2001; Im et al. 2001). Both the de 
Vaucouleurs+exponential disk and Sersic bulge+exponential disk correspond well to visual morphologies, albeit with large scatter; the trend 
between B/T and Hubble type for both profiles is 99% significant (> 2cr) using the Spearman rank test. 
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Fig. 2. — Comparison of Hubble types from FFvDOO to the best fit bulge power n (0.2 < n < 6) for the Sersic bulge+exponential disk 
model; here we consider only members brighter than our cut-off magnitude (mgi4 < 21) with {B/T)c;„ > 0.25. In our cluster sample, we 
find no correlation between Sersic index n and visual type. 
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Fig. 3. — Thumbnail images of 168 cluster members in our sample of 173 fitted galaxies. These images are from the F814W filter Vb)- The 
pixel area in each thumbnail is 15 times the galaxy's isophotal area as defined by SExtractor (typically 10 — 15" on a side). Each set of three 
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F814W magnitude (left), bulge/total luminosity (B/T; middle), and the asymmetry parameter and total fraction of residual light {Ra,Rt', 
right). Higher resolution images of the 168 members are available at http://www.exp-astro.phys.ethz.ch/vy/astronomy 
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Fig. 4. — Distribution of bulge/total luminosities [{B/T)dev] measured with GIM2D in the F814W filter Vq) for 155 cluster members 
("1814 < 21 mags). The dotted line denotes the adopted break between bulge (B/T > 0.4) and disk (B/T < 0.4) dominated galaxies. We 
also determine the {B/T)cicV distribution for the same galaxies in the F606W filter (~ Bo) and find, using the K-S test, that the distribution 
is not dependent on filter. 




i:(Mpc2) 

Fig. 5. — Number fraction of bulge [{B /T)^^y > 0.4; open circles] and disk-dominated [{B /T)^^y < 0-4; open triangles] galaxies as a 
function of local density for the 155 cluster members {m^n < 21 mags); Poisson errorbars are shown. As is observed in lower clusters (Dressier 
1980), the fraction of disk-dominated galaxies increases with decreasing local density. 
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Fig. 6. — Clustering, as defined by the Dressler-Schectman test (Dressier & Shectman 1988), as a function of (B/T)ci^Y for cluster members 
(m8i4 < 21 mags). The scale is 50" ~ 225h~^ kpc. The circles are proportional to how much the local mean redshift and velocity dispersion 
deviate from that of the cluster's; many large circles grouped together indicate significant substructure. As we move from bulge to disk- 
dominated systems (clockwise from upper left panel), we find that bulge-dominated systems are clustered more strongly than disk-dominated 
ones. 
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Fig. 7. — Color-magnitude diagram for different galaxy types as defined by (B /T)deV ■ Open circles are E-SO's (—5 < T < —1), open 
triangles SO-a's (0 < T < 1), open squares spirals and irregulars (2 < T < 15), plus symbols non- typed (mgi4 > 22) members, mergers open 
diamonds (T = 99), and large open circles E+A's. The left panels show {B — V)z vs. B^ (from vD98) for four B/T bins with the most 
bulge-dominated sample {B/T > 0.75) at the top. The CM relation, as determined by vD98 from the E-SO members, is shown as a dotted 
line in each of these panels. The right panels show the corresponding scatter about the CM relation for the four B/T bins. Here C is the 
average deviation from the CM relation and a is the width of the deviation distribution. 
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Fig. 8. — Examples of CL 1358+62 members with a high degree of galaxy asymmetry (Ra > 0.05) and/or high total residual (Rt > 0.1). 
H525 is a galaxy with low R^ and high Rx, while H1475 has high R^ and low R^. The other four galaxies have both R^ > 0.05 and 
Rt ^ 0.1. Higher resolution image is available at http://www.exp-astro.phys.ethz.ch/vy/astronomy/ 
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Fig. 9. — Top Panel: Distribution of {B/T)aeV vs. [OII]A3727 equivalent width {F98) of cluster members {msi4 < 21); negative [Oil] 
EW values correspond to emission. The vertical and horizontal dotted lines refer to B/T = 0.4 and [Oil] EW= — SA, the adopted division 
between bulge/disk-dominated and star-forming ([OII]< — 5A) respectively. Symbols are as in Fig. 7. Bottom Panel: The fraction of galaxies 
in each B/T bin with [Oil] EW< —5 where the dotted line denotes B/T = 0.4. Pure disk systems have a high probability of having significant 
[OII]A3727 emission. 
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Fig. 10. — Top Panels: Correlation between high galaxy asymmetry {Ra > 0.05) and strong [Oil] emission (< — SA) for bulge-dominated 
(top left) and disk-dominated (top right) members; only galaxies with mgi4 < 21 mags are included, and symbols are as in Fig. 7. Vertical 
and horizontal dotted lines show the adopted cuts for high asymmetry and star-forming ([OII]< — 5A) galaxies respectively. The trend 
of increasing [OII]A3727 emission with for the disk- dominated galaxies has 99% confidence with the Spearman rank test. There is no 
detectable trend for the bulge-dominated galaxies. Bottom Panels: Distribution of colors for the two populations where members with no 
significant [Oil] emission are shown as dots, and members with [OII]< — 5A are shown as symbols. Virtually all of the SO-a members with 
significant [Oil] are red. 



20 




Ellipticity cos(i) 

Fig. 11. — Left: The bulge ellipticity of visually typed ellipticals (top) and SO's (bottom; FFvDOO) versus measured (B /T)deV\ here we 
apply an upper magnitude cut so that 19 < m8i4 < 21. We exclude galaxies in the hatched regions [{B/T)^^y < 0.4] as the random error 
associated with ellipticity measurements for these objects is S{Ell) > 0.1. Right: Disk inclination versus {B/T)^^y where only B/T < 0.6 
members are included in the analysis. If E's and SO's truly are two different populations, they should both span the same range in ellipticity 
and inclination. They do not (99% confidence with the K-S test), and there is a conspicuous lack of round {Ell < 0.2) SO's. These points 
strongly suggest these ellipticals are likely to be face-on SO's. 
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Fig. 12. — Top Panel: Difference between input and measured magnitude (Am8i4 = IN — OUT) as a function of measured B/T for a de 
Vaucouleurs bulge+exponential disk in the F814W filter (3600s integration). The sample is split into four apparent (measured) magnitude 
bins with the brightest galaxies in the left-most subpanel. Median values (open circles) and asymmetric errorbars corresponding to la in the 
Am8i4 distributions are shown. Bottom Panel: Same as in the top panel except now compared to the measured half-light radius {ri/2)ouT- 
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Fig. 13. — Comparison of measured versus input values for bulge/total luminosity for a de Vaucouleurs bulge+exponential disk from the 
artificial galaxy catalog in the F814W filter {3600s integration). The y-axis shows A{B/T) = {B/T)ipf - (B/T)ouT- Top Panel: The x-axis 
is measured F814W magnitude, and the sample is split into four {B/T)ouT bins with the most disk-dominated systems in the left-most 
subpanel. Middle Panel: The x-axis is measured B/T value, and the sample is split into four bins with the brightest galaxies in the left-most 
subpanel. Bottom Panel: The x-axis is (fi/2)o(7Ti cut at ^1/2= Sh~^ Kpc to improve the statistics, and the sample is split into four 
bins where the brightest galaxies are in the left-most subpanel. In all three panels, median values (open circles) and asymmetric errorbars 
corresponding to Ic in the A{B/T) distributions are shown. At magnitudes brighter than our cut of mgi4 < 21, systematic errors are 
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Fig. 14. — Left: Histogram showing [A{B /T)deV = {B/T)iiq — {B/T)out] versus measured {B/T)deV for artificial galaxies brighter than 
mgi4 = 21. Right: Comparison of input to measured (B/T)^^y for the same galaxies. 
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[{IN — OUT)/OUT] over the range in {B /T)ouT for ''di a^nd (top, middle, and bottom panels respectively). Asymmetric errorbars 
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radius is robust for {'m,si4)ouT ^ 21 galaxies with systematic and random errors of ~ 10% and ~ 20% respectively. The systematic 
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Fig. 16. — Distribution of input (open symbols) and measured (filled symbols) values for half-light radii (top), disk (middle) and bulge 
(bottom panel) scale lengths as a function of measured {B/T)^^v (bin size 0.1). We include only artificial galaxies with {'msi4)ouT < 21 mags, 
and Icr errorbars are shown. The dotted line denotes the adopted cut between bulge and disk-dominated systems. Recovery of for 
{B/T)ouT < 0.8 galaxies is excellent. Recovery of re for disk-dominated systems is not as good, however, as due to the long wings of the de 
Vaucouleurs profile, flux at large radii can be buried in the background noise. 
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Fig. 17. — Top Panel: The difference between input and measured bulge ellipticity (IN — OUT) for the artificial galaxies (de Vaucouleurs 
bulge+exponential disk). Recovery of bulge elhpticity for these objects is excellent, except for the most disk-dominated systems. The average 
difference is zero, and the associated random la error is A(eZZ) < 0.1. Bottom Panel: The difference between input and measured disk 
inclination for the artificial galaxies. Here, a face-on disk has cos i = I. The recovery of disk inclination for these galaxies also is robust with 
an systematic difference of zero, except for the most bulge-dominated systems; the associated random error is A(cos i) < 0.1. 
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Fig. 18. — Top Panels: The difference between input versus measured value {IN — OUT) of the bulge power n (0.2 < n < 6) from the 
artificial Sersic catalog of {B/T)jf^ > 0.5 galaxies; only galaxies with {'miii4)ouT ^ 21 mags are included. Although 70% of the galaxies 
have differences < 0.5, An can range up to ~ 5 for some. Bottom Panel: The difference in bulge fraction for the same galaxies. Here, the 
distribution is skewed towards underestimating the true bulge fraction where the median difference is ~ 0.14. Many galaxies 25%) now 
are considered disk-dominated systems [{B /T)ouT ^ 0.4] although none have {B/T)j^ < 0.5. 
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Fig. 19. — Left: Histogram showing the difference between input versus measured value (IN — OUT) of B/T for Sersic bulges fitted using 
the r^/^ profile; all of the artificial galaxies have {B/T)jj^ > 0.5. Here we only consider galaxies with measured mgi4 brighter than 21. Right: 
Comparison of input to measured B/T for the same galaxies. As in the case of the Sersic bulge fitted with the Sersic profile, the true bulge 
fraction tends to be underestimated when modeling a Sersic bulge with an r^/* profile. 




Fig. 20. — Left: Histogram showing the difference between input versus measured value {IN — OUT) of B/T for the artificial double 
exponential catalog; only galaxies with (ni8i4)o;7T l£ 21 mags are included. Right: Comparison of input to measured B/T for the same 
galaxies. The obvious problem with using a double exponential profile is that the bulge/disk components arc often reversed, i.e. bulge- 
dominated systems become disk-dominated ones and vice versa. 
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0.61 


21.90 


0.90 


0.49 


0.27 


0.50 


0.01 


0.06 


410 


-2 


0.58 


21.43 


0.87 


0.53 


0.24 


0.64 


0.02 


0.07 


411 





0.50 


21.80 


0.80 


0.58 


0.38 


0.44 


0.01 


0.07 


416 


-2 


0.64 


22.33 


0.84 


0.20 


0.12 


0.20 


0.02 


0.04 


420 


-99 


0.45 


22.71 


0.62 


0.52 


0.35 


0.37 


-0.01 


0.04 
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ID 
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0.46 
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0.55 


1.17 
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0.55 


0.08 


0.20 


594 
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0.60 


22.64 


0.71 


0.13 


0.13 


0.08 


0.02 


0.05 


626 


1 


0.25 


20.94 


0.88 


1.03 


0.29 


0.76 


0.01 


0.07 


1414 


3 


0.65 


22.79 


0.78 


0.45 


0.75 


0.16 


0.01 


0.11 


1455 


99 


0.73 


22.78 


0.82 


0.33 


0.27 


0.26 


0.02 


0.08 


1461 


-99 


0.40 


23.00 


0.77 


0.67 


0.59 


0.41 


0.01 


0.07 


1475 


-2 


0.47 


21.82 


0.84 


0.61 


0.34 


0.48 


0.06 


0.09 


1487 


-99 


0.00 


24.06 


0.46 


0.61 




0.36 


0.03 


0.01 


1524 


-2 


0.16 


22.72 


0.84 


0.43 


0.15 


0.28 


0.01 


0.03 
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Table 1 — Continued 



ID 


Type'' 


IB/T)^ 

\ 1 / 




(B - V)z'' 


\ 1/2/ 


ir'lf 

\ e / 


\ a / 






1563 


-5 


0.79 


22.34 


0.67 


0.62 


0.52 


0.56 


0.03 


0.08 




i 


U.Do 


Z4.00 


U.oU 


U.DD 


U.oi 


u.yo 


n 1 
— u. iz 


n 1 

U.lz 


1616 


-99 


0.29 


22.81 


0.81 


0.46 


0.07 


0.39 


0.01 


0.05 


1630 


-2 


0.00 


22.70 


0.57 


0.46 




0.27 


0.03 


0.09 


1649 


-4 


0.63 


22.37 


0.82 


0.34 


0.24 


0.29 


0.02 


0.07 


1775 


-4 


0.36 


22.38 


0.80 


0.83 


0.58 


0.55 


0.04 


0.12 


1806 


-2 


0.62 


22.59 


0.83 


0.54 


0.55 


0.32 


0.01 


0.08 


1816 


99 


0.69 


23.03 


0.80 


0.48 


0.47 


0.29 


0.01 


0.11 


1829 


-99 


0.12 


24.40 


0.43 


0.44 


0.67 


0.26 


0.05 


0.18 


1842 


-4 


0.11 


23.42 


0.73 


1.06 


0.93 


0.64 


0.01 


0.18 


1865 


-5 


0.69 


22.22 


0.80 


0.54 


0.43 


0.43 


0.01 


0.11 


1871 


5 


0.71 


21.02 


0.89 


0.77 


0.47 


1.04 


0.04 


0.07 


1897 


999 


0.31 


22.34 


0.87 


0.76 


0.67 


0.47 


0.01 


0.11 


1978 


-99 


0.15 


23.23 


0.90 


0.48 


0.24 


0.31 


0.03 


0.10 



'^Hubble types from FFvDOO. Here, morphological types are represented as E (- 
5), SO (-2), Sa (1), Sb (3), Sc(5), Sd (7), Sm (9), and Im (10); mergers are indicated 
with (99), conflicted types as 999, and non-typed galaxies [msiA > 22) with (-99). 

''Measured in the F814W filter using a de Vaucouleurs bulge with exponential 
disk profile. 

'^Converted to rest-frame B and V using transforms from vD98. 



Table 2 

Average [OII]A3727 A Equivalent Widths 



B/T 


N 


[OJ/]A3727 EQ (A) 


771814 


B/T < 0.25 


26 


-11.7± 1.1 


18.3 < m < 20.9 


0.25 < B/T < 0.50 


37 


-1.6 ±0.9 


17.4 <m< 21.0 


0.5 < B/T < 0.75 


60 


-0.3 ±0.8 


17.6 < m < 21.0 


B/T > 0.75 


32 


0.1 ±0.6 


17.8 < TO < 21.0 



